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We calculate the transverse- momentum (pt) distribution for the inclusive hadroproduction of 
B mesons at intermediate values of pT at next-to-leading order (NLO) in a dedicated finite-mass 
scheme using realistic non-perturbative fragmentation functions that are obtained through a global 
fit to e^e~ data from CERN LEPl and SLAG SLC exploiting their universality and scaling vio- 
lations. We find that finite-mass effects moderately enhance the cross section, by about 20% at 



■ Pt = 2mfe, and rapidly fade out with increasing value of pT, so that the zero-mass prediction is 

reached. We also perform comparisons with recent pp data taken by the CDF Collaboration in 
run II at the Fermilab Tevatron and comment on the usefulness of the fixed-flavor-number scheme. 
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I. INTRODUCTION 



Recently there has been much interest in the study of S-meson production in pp colhsions 
at hadron colhders, both experimentally and theoretically. The CDF Collaboration mea- 
sured differential cross sections da/dpx for the inclusive production of B mesons (and their 
anti-particles) in pp collisions at the Fermilab Tevatron as a function of the transverse mo- 
mentum pt in the central rapidity [y) region jl, 2, 3, 4|. The data reported in Ref. [1] were 
collected in the run period from 1992 to 1995 (runs lA and I) at a center-of-mass (cm.) en- 
ergy of \fS = 1.8 TeV and were obtained using fully reconstructed mesons decaying into 
the exclusive final state J/ipK^. The data presented in Ref. come from measurements 
in run II with = 1.96 TeV, where the inclusive differential production cross section of 
J/ ip mesons was used and the fraction of events from the decay of long-lived b hadrons was 
separated by analyzing the lifetime distribution. These b hadrons include B^, B^ , B^, and 
B^ mesons. The data in Ref. were also taken at \/S = 1.96 TeV in run II. In this case, 
the inclusive cross section for the production of B^ mesons was obtained, as in Ref. l|, by 
reconstructing B^ J/ipK^ decays. Very recently, CDF presented preliminary data from 
run II based on events with B^ D^ji^V^ followed by K^n^ and B^ D*^ fi^V^ 

followed by D*^ D^n^ and K^n^ collected with the lepton-plus-displaced-track 

trigger ^. These data explore the range 25 GeV < pj- < 40 GeV for the first time. Although 
the measurement of B mesons is experimentally well defined, theoretical predictions did not 
agree with the data in the past. 

In order to calculate the 5-meson production cross section, the non-perturbative frag- 
mentation function (FF) for the transition b ^ B must be known beforehand. The QCD- 
improved parton model implemented in the modified minimal-subtraction (MS) renormaliza- 
tion and factorization scheme then provides a rigorous theoretical framework for a coherent 
global data analysis. 

In this framework, two distinct approaches for next-to-leading-order (NLO) calculations 
in perturbative QCD have been used for comparisons with experimental data. In the so- 
called massless scheme or zero-mass variable-flavor-number scheme (ZM-VFNS) 
which is the conventional parton model approach, the zero-mass-parton approximation is 
applied also to the b quark, although its mass m is certainly much larger than the asymptotic 
scale parameter Aqcd- In this approach, the b quark is also treated as an incoming parton 
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originating from the (anti)proton, leading to additional contributions besides those from u, 
d, s, and c quarks and the gluon (g). Although this approach can be used as soon as the 
factorization scales associated with the initial- and final-state singularities are above the 
starting scale of the parton distribution functions (PDFs) and the FFs, the predictions are 
reliable only in the region of large pt values, with px ^ m, where terms of the order of 
'm? IPt can safely be neglected. A NLO calculation in this scheme automatically resums 
leading and next-to-leading logarithms (NLL), i.e. terms of the form [a^ ln(p|,/m^)]" and 
ln(pf,/m2)]" with n = 1, 2, 3, . . .. At the same time, all non- logarithmic terms through 
0{as) relative to the Born approximation are retained for m = 351]. 

The other calculational scheme is the so-called massive scheme or fixed-flavor-number 
scheme (FFNS) ^], in which the number of active flavors in the initial state is limited to 
77./ = 4, and the h quark appears only in the final state. In this case, the h quark is always 
treated as a heavy particle, not as a parton. The actual mass parameter m is explicitly taken 
into account along with pt- In this scheme, m acts as a cutoff for the initial- and final-state 
collinear singularities and sets the scale for the perturbative calculations. A factorization of 
these would-be initial- and final-state collinear singularities is not necessary, neither is the 
introduction of a FF for the transition b —>■ B. However, at NLO, terms proportional to 
as ln(pf./m^), where ag is the strong-coupling constant, arise from collinear gluon emissions 
by b quarks or from branchings of gluons into collinear bb pairs. These terms are of order 
0(1) for large values of pt, and with the choice = 0{pt) for the renormalization scale 
they spoil the convergence of the perturbation series. The FFNS with n/ = 4 should thus 
be limited to a rather small range of pr, from Pt = to pt ^ The advantage of this 
scheme is that the m? /p^ power terms are fully taken into account. 

The ZM-VFNS and FFNS are valid in complementary regions of pt-, and it is desirable 
to combine them in a unified approach that incorporates the virtues of both schemes, i. e. 
to resum the large logarithms, retain the full finite-m effects, and preserve the universality 
of the FFs. This is necessary for a reliable and meaningful interpretation of the CDF data 
, which mostly lie in the transition region of the two schemes. An earlier approach 
to implement such an interpolation is the so-called fixed-order-next-to-leading-logarithm 
(FONLL) scheme, in which the conventional cross section in the FFNS is linearly combined 
with a suitably modified cross section in the ZM-VFNS with perturbative FFs, using a 
PT-dependent weight function [ol, [l^. Then the FONLL cross section is convoluted with a 
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non-perturbative FF for the h B transition. These FFs are adjusted to e^e~ data, using 
the same approach, and good agreement with the CDF data was obtained. 

In this work, we wish to present the results of an approach that is much closer in spirit 
to the ZM-VFNS, but keeps all m^/p% power terms in the hard-scattering cross sections. 
This scheme is called general-mass variable-flavor-number scheme (GM-VFNS) and has re- 



cently been worked out for the photoproduction [ll|, [l2| and hadroproduction [3, Q, Q 
of charmed hadrons. In this approach, one starts from the region pt ^ itl and absorbs the 
large logarithms ln(/i|,/m^), where fip is the factorization scale of the initial or final state, 
into the 6-quark PDF of the incoming hadrons and the FF for the 6 — > 5 transition. After 
factorizing the Inm^ terms, the cross section is infrared safe in the limit m — 0, and nj = 5 
is taken in the strong-coupling constant and the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equations. The remaining m-dependent contributions, i.e. the im? /p^ 
power terms, are retained in the hard-scattering cross sections. These terms are very im- 
portant in the region of intermediate pt values, pt ^ fn, and are expected to improve the 
theoretical predictions as compared to the ZM-VFNS. The large logarithms are absorbed 
into the PDFs and FFs by subtraction of the collinearly (mass) singular terms at the initial- 
and final-state factorization scales, respectively. 

It is well known that the subtraction of just the collinearly, i.e. mass singular terms, does 
not define a unique factorization prescription. Also finite terms must be specified. In the 
conventional ZM-VFNS calculation, one puts m = from the beginning, and the collinearly 
divergent terms are defined with the help of dimensional regularization. This fixes the finite 
terms in a specific way, and their form is inherent to the chosen regularization procedure. 
If one starts with m 7^ and performs the limit m — afterwards, the finite terms are 
different. These terms have to be removed by subtraction together with the Inm^ terms in 
such a way that, in the limit — 00, the known massless MS expressions are recovered. 
This matching procedure is needed, since we use PDFs and FFs defined in the ZM-VFNS. 
A subtraction scheme defined in this way is the correct extension of the conventional ZM- 
VFNS to include 6-quark (or similarly c-quark) mass effects in a consistent way. We actually 
include the c-quark contribution in the massless approximation, i.e. we treat the c quark as 
one of the light partons. 

The results of our earlier work on charmed-hadron inclusive production by pp scattering at 
NLO in the GM-VFNS 3,Q, \i\ directly carry over to h hadrons. Then, the h quark is the 
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heavy one, with mass m, while the c quark belongs to the group of light quarks, collectively 
denoted hj q = u, d, s, c in the following, whose mass is put to zero. Furthermore, we need 
PDFs and FFs implemented with = 5 in the MS factorization scheme. Non-perturbative 
FFs for the transitions a B^, where a = g,q, q, b, b, were extracted at leading order (LO) 
and NLO already several years ago jsj using data for the scaled-energy (x) distribution da/dx 
of e+e" ^ S + X at = 91.2 GeV measured by the OPAL Collaboration at CERN LEPl 
Q. 

We note that our implementation of the GM-VFNS is similar to the Aivazis-Collins- 



Olness-Tung (ACOT) [17|] scheme formulated for the initial state of fully inclusive deep- 
inelastic scattering. The extension of this scheme to the inclusive production of heavy 
partons was considered in Ref. |1S], where the resummation of the final-state collinear loga- 
rithms was only performed to LO and parton-to-hadron FFs were not included. A discussion 
of the differences between our approach and the one in Ref. tl8| concerning the collinear 

n 

subtraction terms can be found in Ref. |14j |. 

This paper is organized as follows. In Sec. [Tll we introduce new NLO sets of 5-meson FFs. 
In Sec. mil we numerically analyze the GM-VFNS predictions with regard to the impact of 
the m-dependent terms and the relative importance of the various partonic initial states. In 
Sec. llVt we compare the predictions of the GM-VFNS, and also those of the ZM-VFNS and 
FFNS, with CDF data from run II . Our conclusions are contained in Sec. |Vl 

II. NON-PERTURBATIVE B-MESON FRAGMENTATION FUNCTION 

As input for the calculation of inclusive i?-meson production cross sections one needs a 
realistic non-perturbative FF describing the transition of the b (b) quark into a B meson. 
Such a FF can be obtained only from experiment. In Ref. [5^ LEPl data for the distribution 
in the scaled 5-meson energy, x = 2Eb/VS, from OPAL [16] were fitted at LO and NLO in 
the ZM-VFNS using three different ansaetze for the b B FF at the starting scale fip = fJ-o 
of the DGLAP evolution, including the ansatz by Peterson et al. (lol, 

and the simple power ansatz j2o| |. 

D{x,^ll)=Nx''{l-xY. (2) 



The best fit was obtained for tlie Peterson ansatz. In Ref. j^, tlie starting scale was taken 
to be fiQ = 2m, witli m = 5.0 GeV. Tlie a ^ B FFs for a = g,q,q were assumed to be zero 
at fip = f^Q and generated through the DGLAP evolution to larger values of fip. 

In the meantime, new and more precise measurements of the cross section of inclusive 
5-meson production in e^e~ annihilation on the Z-boson resonance have been published 
by the ALEPH ^21|, OPAL 22], and SLD [9] collaborations, which motivates us to update 
the analysis of Ref. This also gives us the opportunity to adjust some of the choices 
made in Ref. sl, to conform with the conventions underlying modern PDF sets. In fact, 
for our numerical analysis, we use the NLO proton PDF set CTEQ6.1M, based on the MS 
prescription, by the Coordinated Theoretical-Experimental Project on QCD (CTEQ) 24 1. 
In this set, the 6-quark PDF has its starting scale ai fio = m with m = 4.5 GeV. The 
mass values used in PDFs and FFs have, of course, to be chosen consistently in order to 
avoid the appearance of terms proportional to ln(/iQ/m^) in the NLO corrections. While a 
shift in the starting scale from fiQ = 2m with m = 5.0 GeV to fio = m with m = 4.5 GeV 
changes the 6-quark FFs only marginally at fip values relevant for the e'^e~ annihilation 
cross sections to be used in the fit, it does have a significant effect on the g ^ B FF, which 
greatly affects the cross section predictions for the Tevatron. The size of the analogous effect 
for D*^ FFs is investigated in Fig. 1 of Ref. 25|, which uses /xq = f^c through comparison 



with Ref. 

sets 



21 



26 L which uses /iq = 2mc- We thus perform a combined fit to these three data 
23( 1 using fiQ = m with m = 4.5 GeV as in Ref. 2J]. Furthermore, we adopt 
from Ref. 24 1 the NLO value A^ = 227 MeV appropriate for = 5, which corresponds to 



2l 



a 



(5), 



TTiz) = 0.1181. As in Ref. |5|], we chose the renormalization and factorization scales to 
he = lip = VS- We use the ansaetze of Eqs. ([1]) and ([2]) for the b B FF at /ij? = /iq, 
while the g,q —>■ B FFs are taken to vanish at fip = yUo and are generated through the 
DGLAP evolution. In order to obtain acceptable fits, we have to omit some of the data 
points in the small-x region. Specifically, we only include the ALEPH data with x > 0.25, 
the OPAL data with x > 0.325, and the SLD data with x > 0.28. At the other end of the 
X range, we include all data points up to x = 1. Altogether we use 18, 15, and 18 data 
points of the ALEPH, OPAL, and SLD sets, respectively. Since we only include in the fit 
data from the Z-boson resonance, finite-m effects, being of relative order m^/m| = 0.2%, 
are greatly suppressed, so that we are comfortably within the asymptotic regime where the 
GM-VFNS is equivalent to the ZM-VFNS. 
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TABLE I: Fit parameters of the 6-quark FFs in Eqs. ([T]) and ([2]) at the starting scale fio = m = 
4.5 GeV and values of Xd o f achieved. All other FFs are taken to be zero at /xq = ^■ 

Nap e xVd.o.f 

0.06634 - - 0.008548 21.37 

4684.1 16.87 2.628 - 1.495 



TABLE II: Branching fractions B{fj,p) and average energy fractions {x){fj,F) evaluated at factor- 
ization scales lip = 4.5, 9.0, and 91.2 GeV using the h ^ B FFs based on the Peterson and power 
ansaetze. 



(in GeV) 


B{^XF) 


Peterson 


B{^lF) 


power 

{x){flF) 


4.5 


0.3994 


0.8098 


0.4007 


0.8312 


9.0 


0.3935 


0.7542 


0.3955 


0.7730 


91.2 


0.3767 


0.6403 


0.3803 


0.6537 



The values of the parameters in Eqs. ([T]) and ([2]) obtained through the fits based on the 
Peterson and power ansaetze are listed in Table [T] together with the respective values of 
per degree of freedom, x^/d.o.f. The corresponding da/dx distributions are compared with 



the ALEPH 



2l|, OPAL |22|, and SLD |23| data in Fig. [H These three data sets mostly 



overlap and can hardly be distinguished in the figure. We observe from Table [T] and Fig. [T] 
that the power ansatz yields an excellent overall fit to the selected data points. There are 
deviations at x < 0.3, which are due to the exclusion of data points from the fit. The 
X^/d.o.f value for the combined fit is 1.495. The individual x^/d.o.f values of the ALEPH, 
OPAL, and SLD data sets are 0.861, 2.350, and 1.410, respectively. On the other hand, the 
Peterson ansatz leads to an intolerable description of the data, yielding x^/d.o.f = 21.37 
for the combined fit and similar values for the individual data sets. This ansatz has only 
two free parameters, N and e, and is just not flexible enough to account for the very precise 
experimental data. 

Besides the h ^ B FF itself, also its first two moments are of phenomenological interest 
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FIG. 1: Comparisons of the ALEPH [21] (circles), OPAL [22] (squares), and SLD 
data with the NLO fits using (a) the Peterson ansatz ([1]) and (b) the power ansatz ([2|). The 
initial factorization scale for all partons is /xq = w. = 4.5 GeV. The different symbols can be better 
distinguished in the electronic edition, where the figures can be enlarged. 
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and subject to experimental determination. They correspond to the b ^ B branching 
fraction, 

B{hp) = / dxD{x,jj%), (3) 

and the average energy fraction that the B meson receives from the b quark, 

1 



= / dxxD{x,^l), (4) 

where the cut a;cut = 0.15 excludes the x range where our formalism is not valid. We 
observe from Table [TTl that the Peterson and power ansaetze lead to rather similar results 
for B{hf) and {x){fiF), the ones for the power ansatz being slightly larger. While B{hf) is 
practically independent of {x){fJ'F) is shifted towards smaller values through the evolution 
in fiF- It is interesting to compare the results for (x)(91.2 GeV) in Table HTl with the values 



quoted by ALEPH, OPAL, and SLD,^ which read 0.7361 ± 0.0061 (stat) ± 0.0056 (syst) |2l|, 



0.7193±0.0016 (stat)™* (syst) [2^, and 0.709±0.003 (stat)±0.003 (syst)±0.002 (model) 



23l |. respectively. We observe that the experimental results lie systematically above ours. 
However, one must keep in mind that the experimental results refer to the first moment 
of the measured cross section distribution da/dx, which naturally includes all orders and 
also contributions from gluon and light-quark fragmentation, while ours are evaluated from 
the 6 ^ 5 FF at NLO in the MS scheme via Eq. @). Of course, the 6 ^ 5 FF and its 
moments depend on scheme, order, and implementation issues such as the functional form 
of the ansatz at the starting scale /iq and the value of /xq itself, and thus do not represent 
physical observables by themselves. Nevertheless, comparisons of the quantities B{^f) and 
{x){iif) defined in Eqs. ([3]) and (jl]), respectively, with their experimental counterparts are 
useful to check the dominance oib B fragmentation and are routinely performed in the 
literature (see, e.g., Ref. [^I). 

III. THEORETICAL PREDICTIONS FOR pp ^ B + X 

We are now in a position to perform a numerical analysis. We consider the inclusive 
cross section of pp B + X, where B stands for the average of the B^ and B~ mesons, 
at = 1.96 TeV as in run II at the Tevatron. We concentrate on the px distribution 
integrated over |?/| < 1 cor resp onding to the central region of the CDF detector. We use the 



CTEQ6.1M proton PDFs 



24l | and the 5-meson FFs based on the power ansatz presented in 



Sec. [ni both implemented at NLO with = 227 MeV and m = 4.5 GeV. For simphcity, we 
use a common factorization scale for the initial and final states. We set the renormalization 
and factorization scales to fiR = ^RrriT and = ^Ff^T, where m-r = \/ + m? is the 
transverse mass of the h quark and ^r and ^f are introduced to estimate the theoretical 
uncertainty. Unless otherwise stated, we use the default values = = ^- With our 
default choices fiQ = m and /i_p = ""^t, we have /i^ — > /io as ^ 0. In this limit, the FFs 
and 6-quark PDF should fade out and quench the cross section, leading to a turn-over of the 
Pt distribution. However, the precise location of the maximum and other details of the line 
shape are also subject to other implementation issues of the GM-VFNS. We shall return to 
this topic in Sec. HVl 

The calculation of the cross section (Pa/{dpTdy) of -B-meson hadroproduction at NLO 
in the GM-VFNS proceeds analogously to the case of D mesons outlined in Ref. 13|. Now, 
m denotes the mass of the h quark, and the c quark belongs to the group of light quarks g. 



(a) 



(b) 



B 



g 



n 




(c) 



FIG. 2: Examples of Feynman diagrams leading to contributions of (a) class (i), (b) class (ii), and 
(c) class (iii). 
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whose mass is put to zero. The NLO cross section consists of three classes of contributions. 

1. Class (i) contains all the partonic subprocesses with a b,b ^ B transition in the final 
state that have only light partons {g, q, q) in the initial state, the possible pairings 
being gg, gq, gq, and qq. A Feynman diagram representing this class is shown in 
Fig. da). 

2. Class (ii) contains all the partonic subprocesses with a.b,b ^ B transition in the final 
state that also have 6 or 6 quarks in the initial state, the possible pairings being gb, 
gb, qb, qb, qb, qb, and bb [see Fig. [2](b)]. 

3. Class (iii) contains all the partonic subprocesses with a g,q,q B transition in the 
final state [see Fig. [2](c)]. 

In the FFNS, only the contribution of class (i) is included, but the full m dependence is 



retained 



Isj. On 



the other hand, in the ZM-VFNS, the contributions of all three classes are 



taken into account, but they are evaluated for m = [27j. In the GM-VFNS, the class- (i) 
contribution of the FFNS is matched to the MS scheme, through appropriate subtractions 
of would-be collinear singularities, and is then combined with the class- (ii) and class- (iii) 
contributions of the ZM-VFNS; thus, only the hard-scattering cross sections of class (i) 
carry explicit m dependence. Specifically, the subtractions affect initial states involving 
g ^ bb splittings and final states involving g ^ bb, b ^ gb, and b ^ gb splittings, and 
they introduce logarithmic dependences on the initial- and final-state factorization scales 
in the hard-scattering cross sections of class (i), which are compensated through NLO by 
the respective factorization scale dependences of the 6-quark PDF and the b ^ B FF, 
respectively. The explicit form of the subtractions may be found in Ref. p^]. A certain 
part of the class- (ii) and class- (iii) contributions is due to Feynman diagrams with internal 
6-quark lines; another one is due to diagrams with external 6-quark lines and contains m- 
dependent logarithms, which are resummed. In the FFNS, the m dependence of these 
contributions would only enter beyond NLO, which is reflected in the ZM-VFNS by the 
generic suppression of the 6-quark PDF relative to the gluon and g-quark ones and of the 
gluon and g-quark FFs relative to the 6-quark one. This entitles us to omit this m dependence 
by calculating the contributions of classes (ii) and (iii) in the ZM-VFNS. It turns out that q- 
quark fragmentation contributes negligibly. However, the gluon fragmentation contribution 
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reaches approximately 50% at small values oipr, and its relative contribution decreases only 
rather mildly towards larger values of pr- 

We first investigate the effect of the finite-m terms in the hard-scattering cross sections 



10^ E 



10^ r 



102 r 



10 ^ 



1 ^ 



10-1 



m = 



m = 4.5 GeV Xv 



10 



LO 

gg ^ bb -\- qq 



1 



15 



bb - 



20 25 
PT[GeV] 



FIG. 3: Transverse-momentum distribution da/dpT oipp^ B + X cm. energy = 1.96 TeV 
integrated over the rapidity range \y\ < 1. The contributions of class (i) evaluated at LO in the 
ZM-VFNS (dashed line) and the GM-VFNS (solid line), but with the NLO versions of a^, the 
PDFs, and the FFs, are compared. 
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and thus concentrate on the contribution of class (i) for the time being. This effect can 
already be studied at LO, where the partonic subprocesses read g+g b+b and q+q b+b. 
To this end, we simply switch off the NLO terms in the hard-scattering cross sections while 
keeping a^, the PDFs, and the FFs at NLO, although, strictly speaking, this does not 
represent a genuine LO analysis. The results are shown in Fig. 121 where the dashed and 
solid lines refer to the results for zero and finite values of m, respectively. We observe that 
these results rapidly approach each other with increasing value of pt- At pr = 7.5 GeV, 
the finite-m result is 33% smaller than the m = one, a relative difference of the order of 
rn? Ip^, as expected. 

We now turn to NLO by switching on the QCD corrections to the hard-scattering cross 
sections of class (i). The results for m = and finite m are shown in Fig. [Has the upper 
and lower solid lines, respectively. They constitute parts of the final ZM-VFNS and GM- 
VFNS results. In both cases, the contributions of classes (ii) and (iii) for m = still 
must be added to obtain the full predictions to be compared with experimental data. The 
class- (i) contributions in the ZM-VFNS and GM-VFNS schemes are, therefore, entitled to 
be negative and they indeed are, for px ^ 76 GeV and pt ^ 10 GeV, respectively, as 
may be seen from Fig. |H Comparing the ZM-VFNS and GM-VFNS results, we notice 
that the finite-m effects are significant for px ^ 10 GeV and even cause a sign change for 
10 GeV < Pt ^ 76 GeV. However, as will become apparent below, the contributions of class 
(i) are overwhelmed by those of classes (ii) and (iii), so that the finite-m effects are washed 
out in the final predictions, except for very small values of px- It is instructive to study 
the relative importance of the (/(/-initiated contributions. They are also included in Fig. H] 
for m = and finite m as the upper and lower dashed lines, respectively. They exhibit a 
similar pattern as the full class- (i) contributions and dominate the latter in the small-pr 
range. Comparing Fig. Hlwith Fig. 2(c) in Ref. 13(], we observe that the relative influence 
of the finite-m effects is much smaller in the c-quark expected because the c quark 

is much lighter than the b quark. 

In the remainder of this section, we work in the GM-VFNS and also include the contri- 
butions from classes (ii) and (iii), i.e. we allow for b (anti)quarks in the initial state and 
diQiQ ^ B fragmentation. It is interesting to study the relative importance of the various 
initial states. In Fig. [5], the total result in the GM-VFNS (solid line) is broken up into the 
contributions from initial states consisting of (1) one gluon and one b (anti) quark (upper 
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5 10 15 20 25 

PT[GeV] 

FIG. 4: Transverse-momentum distribution da/dpT of — ^ S + X at cm. energy \/S = 1.96 TeV 
integrated over the rapidity range \y\ < 1. The contributions of class (i) (sohd hues) and their gg- 
initiatcd parts (dashed lines) evaluated at NLO in the ZM-VFNS (upper lines) and the GM-VFNS 
(lower lines) are compared. 

dashed line); (2) one q (anti)quark and one h (anti)quark (middle dashed line); (3) two h 
(anti)quarks (lower dashed line); (4) one gluon and one q (anti)quark or two q (anti)quarks 
(lower dotted line); and (5) two gluons (upper dotted line). If it were not for the class- (iii) 
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PT[GeV] 

FIG. 5: Transverse-momentum distribution da/dpr oipp ^ B + X at cm. energy ^ = 1.96 TeV 
integrated over the rapidity range \y\ < 1. The total NLO result in the GM-VFNS including classes 
(i)-(iii) (solid line) is broken up into the contributions from initial states consisting of (1) one gluon 
and one b (anti)quark (upper dashed line); (2) one q (anti)quark and one b (anti)quark (middle 
dashed line); (3) two b (anti)quarks (lower dashed line); (4) one gluon and one q (anti)quark or 
two q (anti)quarks; and (5) two gluons. 
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contribution, then the combination of contributions (4) and (5) would coincide with the 
class-(i) contribution, considered in Fig. HI and the combination of contributions (l)-(3) 
would coincide with the class- (ii) contribution. However, in Fig. the class- (iii) contribu- 
tion is distributed among the contributions (l)-(5) according to the respective initial states. 
We observe from Fig. O that the partonic subprocesses with one b or b quark in the initial 
state make up the bulk of the cross section throughout the entire mass range considered. 
Specifically, the contribution from the subprocesses where the second incoming parton is 
a gluon (1) is more than twice as large than the one where this is a light (anti)quark (2), 
and it is even larger than the purely gluon-initiated contribution (5), which is a surprising 
finding in view of the enormous gluon luminosity in pp collisions at a cm. energy of almost 
2 TeV. On the other hand, the contribution from two incoming b (anti) quarks (3) is greatly 
suppressed, being less than 1% of the full result. The contribution due to light-parton initial 
states with no more than one gluon (4) ranks between contributions (2) and (3), and it 
is negative for pr ^ 7 GeV. As explained above, the difference between the ^f^f-initiated 
GM-VFNS contributions in Figs. H] (lower dashed line) and [5] (upper dotted line) is due to 
QiQiQ ^ B fragmentation being included in the latter. Obviously, this additional contribu- 
tion is quite significant throughout the whole px range considered. Comparing the difference 
between the class- (i) contributions in the ZM-VFNS and the GM-VFNS in Fig. Hlwith the 
total result in Fig. O we anticipate that the finite-m effects on the latter will be rather 
moderate, except for very small values of pt-, where the ZM-VFNS is expected to break 
down anyway. Also taking into account that the class- (i) contributions considered in Fig. H] 
are less negative (or even positive) in the GM-VFNS than they are in the ZM-VFNS, we 
conclude that the finite-m effects will moderately enhance the cross section. This point will 
be subject to further investigation in the next section. 



IV. COMPARISON WITH CDF DATA 



We are now ready to compare our NLO predictions for the cross section distribution 
da/dpT with Tevatron data. We focus our attention on the more recent CDF data from 
run II published in Refs. {2, 3]. This comparison is presented for the GM-VFNS in Fig. [6l 
where the solid line represents the central prediction, for C,r = C,f = ^, and the dashed lines 
indicate the maximum and minimum values obtained by independently varying C,r and 
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FIG. 6: Transverse-momentum distribution da/dpT of pp B + X at cm. energy yS* = 1.96 TeV 
in.eg„ted „ve. the ™p.d., ™nge H < 1. The cent.^ NLO p.e«.n w.th & ^ ^1 (.ol.d 

line) of the GM-VFNS is compared with CDF data from Refs. [2] (open squares) and [3] (solid 
squares). The maximum and minimum values obtained by independently varying and in the 
range 1/2 < $,r,£,f ^ 2 with the constraint that 1/2 < ^r/(,f ^ 2 are also indicated (dashed lines). 

in the range 1/2 < C,r,C,f < 2 with the constraint that 1/2 < C,r/C,f < 2. The maximum and 
minimum values correspond to .^i? = 2 and C,f = 1/2? respectively. The variation with C,r is 
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considerably milder than the one with and only leads to a modest broadening of the error 
band. For C,f < ^, f^F reaches the starting scale (10 = 171 for the DGLAP evolution of the FFs 
and the 6-quark PDF at pt = ^a/I/'Cf ~ 1- Fo^' smaller values of pt, there is no prediction 
because the FFs and the 6-quark PDF are put to zero for < A^o- This explains why the 
Pt distribution for = 1/2 only starts at pt = V^m ^ 7.8 GeV. The most recent data ^ 
nicely agree with the GM-VFNS result. In fact, they lie close to the central prediction, with 
a tendency to fall below it in the lower p^ range, and they are comfortably contained within 
the theoretical error band. Obviously, the notorious Tevatron S-meson anomaly, with data- 
to-theory ratios of typically 2-3 |l| , that has been with us for more than a decade has finally 
come to its end, thanks to both experimental and theoretical progress (see also Ref. 28| for 
a recent status report on the observed and predicted cross sections at the Tevatron). The 
previous CDF data [2], based on a measurement of J/ip + X final states, are compatible 
with the latest ones for < 12 GeV, but systematically undershoot them for larger values 
of Pt- This potential inconsistency becomes even more apparent by noticing that Fig. [6] 
only contains 4 out of the 13 data points for pt > 12 GeV quoted in Ref. and that the 
omitted data points neatly line up with the selected ones. This possibly suggests that the 
systematical errors in Ref. [2,], and perhaps also in Ref. [3'], might be underestimated and 
that the overall normalization might need some adjustment. Incidentally, the preliminary 
CDF data {4] fall right in the middle between those from Refs. and Sj. 

The measured pt distribution of Ref. 2|] reaches down to = and exhibits a maximum 
at ~ 2.5 GeV. As we shall see below, this small-pr behavior is correctly reproduced in the 
FFNS without DGLAP-evolved FFs, which only receives contributions of class (i) without 
any subtractions. It is clear that our present implementation of the GM-VFNS is not suitable 
for cross section calculations in the small-p^ region. Although the GM-VFNS is designed 
to approach the FFNS and the ZM-VFNS in its regions of validity without introducing 
additional ad-hoc matching factors, to implement this numerically is a non-trivial task due 
to necessary cancellations between different terms in the calculation. Stable computer codes 
including these features have been developed for the fully inclusive case and are used in 
global analyses of proton PDFs, e.g. in the CTEQ studies [2A] using the ACOT scheme 
For one-particle inclusive processes, the problem to achieve such cancellations is complicated 
by the extra factorization scale; to obtain a smooth transition from the GM-VFNS to the 
FFNS, one has to carefully match terms that are taken into account at fixed order with 
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terms that are resummed to higher orders in the PDFs and FFs. In addition, it remains to 
be investigated whether a proper scale choice in the small-pr range is required and helpful 
to ensure that the FFs and 6-quark PDF are sufficiently suppressed already at pt = 0{m). 

The GM-VFNS prediction in Fig. [H] exhibits a sizeable scale uncertainty for < 2m. 
As mentioned above, the pt distribution for = 1/2 only starts at pt ~ 7.8 GeV. These 
undesirable features will eventually be removed once the matching with the FFNS is specified 



and implemented. This is a so-called implementation issue [29| that needs to be added on top 
of the definition of the pure GM-VFNS. This is beyond the scope of the present paper, which 
is concerned with the intermediate px range, and will be treated in a future publication. At 
this point, we would like to recall how this implementation issue is handled for the FONLL 
scheme jOj. In that scheme, the FFNS and ZM-VFNS calculations are merged in such a 
way that the contribution that is added on to the FFNS result, i.e. the ZM-VFNS result 
with the zero-mass limit of the FFNS result subtracted, is multiplied by a weight function 
of the form Pt/{Pt + c^'^^) with c = 5 to model a smooth transition. Furthermore, the 
variable pt of the subtracted ZM-VFNS contribution is shifted to become m^. In the region 
where the ZM-VFNS prediction has a large scale uncertainty, i.e. where pt is 2-3 times 
larger than m say, this weight function is still rather small, ranging from 14% to 26%. Thus, 
this weight function not only smoothens the transition, but also ensures that the sizeable 
theoretical uncertainty of the ZM-VFNS component in the transition region is not reflected 
in the FONLL prediction, creating the impression that the latter has a small theoretical 
error. Of course, this source of theoretical uncertainty unavoidably resurfaces when the 
form of the weight function, which is a priori unknown, is varied, e.g. by changing the value 
of its parameter c. Unfortunately, such a variation is not included in the theoretical error of 
recent FONLL predictions Q, 34|- 

We now extend our numerical analysis to include the NLO prediction in the FFNS, with 
Hf = 4 massless quark flavors in the initial state, which allows us to also compare with the 
small-pT data from Ref. [J]. In the FFNS analysis, we evaluate a'f'^\j2ji) wit h Uf = 4 and 
= 326 MeV 3], while we continue using the CTEQ6.1M proton PDFs 24|, in want 
of a rigorous FFNS set with nf = 4. In the FFNS, there is no room for DGLAP-evolved 
FFs, and only b,b ^ B transitions are included. For simplicity, we identify b (anti)quarks 
with B mesons and account for non-perturbative effects by including the branching fraction 



B{b ^ B) = 39.8% 30| as an overall normalization factor, i.e. we use a 6 — > 5 FF of the 
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FIG. 7: Transverse-momentum distribution da/dpx of pp B + X at cm. energy V'S' = 1.96 TeV 
integrated over the rapidity range \y\ < 1. The central NLO predictions in the FFNS with Uf = 4 
and without FFs (dot-dashed hue), the ZM-VFNS (dashed hne), and the GM-VFNS (sohd hne) 



are compared with CDF data from Refs. [2] (open squares) and 
the historical FFNS prediction, evaluated with PDF set MRSDO 



Jsi (solid squares). For reference, 
311], ab ^ B FF of Peterson type 



ig[ with e = 0.006, rrib = 4.75 GeV, and = 215 MeV, is also shown. 



form D{x) = B{h B)5{1 — x), while the g,q,q B FFs are put to zero. In Fig. [71 
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the central FFNS (dot-dashed hne), ZM-VFNS (dashed hne), and GM-VFNS (sohd hne) 
predictions, for = = ^, are compared with the CDF data from Refs. As in 

Fig. El some of the data points with px > 7 GeV from Ref. are omitted for clarity. Since 
the ZM-VFNS and our present implementation of the GM-VFNS are not applicable to the 
small-pr range, we show the respective predictions only for > 2m = 9 GeV. The GM- 
VFNS prediction shown in Fig. [7] is identical with the central one in Fig. [61 By construction, 
it merges with the ZM-VFNS prediction with increasing value of pt- In accordance with 
the expectation expressed in the discussion of Figs. [3] and HI the difference between the 
GM-VFNS and ZM-VFNS results is rather modest also at pt ^ 2m, since the m-dependent 
contribution, of class (i), is numerically small and overwhelmed by the m-independent ones, 
of classes (ii) and (iii). The FFNS prediction faithfully describes the peak structure exhibited 
by the next-to-latest CDF data 21] in the small-pr range and it also nicely agrees with the 
latest CDF data [s] way out to the largest pt values. In fact, for px > 4m, where its 
perturbative stability is jeopardized by unresummed logarithms of the form ln(m|./m^) > 3, 
the FFNS prediction almost coincides with the GM-VFNS one, where such large logarithms 
are resummed. This is a pure coincidence, which becomes even more apparent if we also 
recall that the implementation of the b,b ^ B transition in the FFNS is not based on a 
factorization theorem and quite inappropriate for such large values of pt- 

In Figs. [6] and U\ we limited our considerations to the range px < 25 GeV, where the 
pubhshed CDF data B,!^ are located. However, the preliminary CDF data 1^, collected in 
the very central part of the detector {\y\ < 0.6), cover the range 9 GeV < < 40 GeV, and 
it is interesting to confront them with the NLO predictions of the three schemes considered 
here. Moreover, it is instructive to study the breakdown of the FFNS at sufficiently large 
values of pt due to unresummed large logarithms. For these purposes, we show in Fig. [8] 
an extension of Fig. [7] including the preliminary CDF data. We observe that the GM- 
VFNS result steadily merges with the ZM-VFNS one as the value of px is increased, the 
relative deviations being 11%, 6%, and 3% at pr = 20, 30, and 50 GeV, respectively. 
The FFNS result breaks even with the GM-VFNS one at about pt = 20 GeV (see also 
Fig. E]) and exceeds the latter for larger values of pt, the relative deviations being 20%, 
48%, and 100% at pt = 30, 50, and 100 GeV, respectively. Our results indicate that a 
measurement of the pt distribution up to 40-50 GeV at the Tevatron could be able to 
resolve the difference between the FFNS and the VFNS and so to establish for the first time 
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FIG. 8: Transverse-momentum distribution da/dpT olpp ^ B + X at cm. energy ^/S = 1.96 TeV 
integrated over the rapidity range \y\ < 0.6. The central NLO predictions in the FFNS with nj = 4 
and without FFs (dot-dashed line), the ZM-VFNS (dashed line), and the GM-VFNS (solid hne) 
are compared in the large-pr range. 
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the break-down of the FFNS due to unresummed logarithms in the inclusive hadroproduction 
of heavy hadrons. This important question deserves a careful examination of the theoretical 
uncertainties, which we leave for a future publication. The preliminary CDF data point 
in the bin 29 GeV < < 40 GeV favors the ZM-VFNS and GM-VFNS resuhs, while it 
undershoots the FFNS result. 

We conclude this section with an interesting observation that, in retrospect, sheds some 
new light on the Tevatron 5-hadron anomaly mentioned above and does not appear to be 
sufficiently well known to the community. In fact, the common perception that the CDF 
data generally overshoot the FFNS prediction, frequently denoted as NLO QCD 

in the literature, by a factor of 2-3 is entirely due the use of obsolete theoretical input. 



In fact, the FFNS prediction that has 



Deen serving as a benchmark for some 15 years and 
28] is evaluated with the proton PDF set MRSDO 



still does even in very recent papers 

by Martin, Roberts, and Stirling |3l|, which has been revoked by these authors. It has an 
unacceptably weak gluon and a small value of A^, namely = 215 MeV translating 
into a^^\mz) = 0.111, which is 3.3 standard deviations below the present world average 
af'\mz) = 0.1176 ± 0.0020 [3^. Other inputs include = 4.75 GeV and a Peterson 
FF parameter of e = 0.006, extracted from a fit to e+e^ annihilation data from the pre- 
EP/SLC era using a Monte-Carlo event generator based on massless LO matrix elements 



32] • For reference, the historical FFNS prediction evaluated with this choice of input is 
also included in Fig. [7l Since the FONLL prediction [l^ is designed to merge with the 
FFNS one at low values of px, the additional contribution being faded out by a weight 
function of the form Px/{Pt + 25m^), the striking gap between the historical NLO QCD 
prediction and the FONLL prediction, based on up-to-date input information, in Fig. 10 of 
Ref. {3}] impressively illustrates the advancement in the PDF and as determinations. The 
tuning of FFs in connection with the resummation of leading and next-to-leading logarithms, 
emphasized in the second paper of Ref. ^|, is actually of minor importance. 



V. CONCLUSIONS 



For several years, the 5-meson production rates measured at DESY HERA, CERN LEP2, 
and Tevatron have been notoriously exceeding, by up to a factor of three, the usual NLO 
QCD predictions for massive b quarks, i.e. those in the FFNS (5-meson anomaly). This 
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has even triggered theoretical attempts to interpret this deviation as a signal of new physics 
beyond the standard model jssl. However, it remained to be clarified if this deviation 
could be explained by improving and refining the QCD prediction itself. In this connection, 
two of us, together with Binnewies, pointed out almost a decade ago that the ZM-VFNS 
provides a rigorous theoretical framework for a coherent study of S-meson production in 
high-energy e^e~, pp, and other collisions, since the factorization theorem guarantees the 
universality of the i?-meson FFs In fact, the ZM-VFNS prediction js], Ql was found to 
nicely agree with the CDF data from Tevatron runs lA and I [1]. However, a necessary 
condition for the applicability of the ZM-VFNS is that the energy scale that separates 
perturbative hard scattering from non-perturbative fragmentation (final-state factorization 
scale fip) is sufficiently large compared to the 6-quark mass m, and it had never been quite 
clear how large the ratio fip/rn actually needed to be in order for finite-m effects to be 
negligible. In fact, the authors of Ref. [3j| asserted in a footnote that mass corrections have 
a large size up to pr ~ 20 GeV and that "lack of mass effects [sl will therefore erroneously 
overestimate the production rate at small pt" 

In the present paper, we addressed this problem by performing a comparative analysis of 



i?-meson hadroproduction in the ZM-VFNS and the GM-VFIN 
applied to D-meson production in 77 [11], ep |l2], and pp 13, 



S, which we had successfully 



For this, we also updated the determination of 5-meson FFs p] by fitting to recent e^e 



data from ALEPH 



2l|, OPAL 



and SLD 



1141 



15| collisions in the past. 



23| and also adjusting the values of m and the 



energy scale fio where the DGLAP evolution starts to conform with modern PDF sets 24 1. 
We found that finite-m effects moderately enhance the pt distribution; the enhancement 
amounts to about 20% at px = 2m and rapidly decreases with increasing value of pt, falling 
below 10% at pt = 4:771. This finding contradicts earlier claims 3J] in all respects. Such 
effects are comparable in size to the theoretical uncertainty due to the freedom of choice in 
the setting of the renormalization and factorization scales. For comparison, we also evaluated 
the Pt distribution in the FFNS, with n/ = 4, using a delta-function-type h ^ B FF without 
DGLAP evolution. 

^ n 

Confronting the three NLO predictions with the latest [3] and next-to-latest |^ CDF data 
sets published, we found that all of them agree rather well with the latest one, with pt > 
7 GeV. Despite unresummed large logarithms and poorly implemented fragmentation, the 
FFNS prediction happens to almost coincide with the GM-VFNS one in the range 15 GeV < 



24 



Pt ^ 25 GeV. The FFNS prediction also nicely reproduces the peak exhibited about pt ~ 
2.5 GeV by the next-to-latest CDF data [2]. By contrast, the historical benchmark result 
based on obsolete proton PDFs and a value of af'\mz) falling short of the present world 
average by 3.3 standard deviations, which goes under the name NLO QCD in the literature 
and is used as a reference point even in most recent papers 3|, |28|], undershoots the CDF data 
by the familiar factor of 2-3. This illustrates that the progress in our understanding of the 
proton PDFs and our knowledge of ai^\mz) is instrumental in overcoming the long-standing 
Tevatron 5-hadron anomaly in the low to intermediate px range. The preliminary CDF data 
4 favor the ZM-VFNS and GM-VFNS results in the upmost bin, 29 GeV < pt < 40 GeV, 
while they undershoot the FFNS result. 

It is desirable to extend the applicability of the GM-VFNS down to pt = 0- This requires 
matching with the FFNS. To achieve this in a way that avoids ad-hoc weight functions is a 
non-trivial task and is left for future work. 
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